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A  RUNOFF  HYDROGRAPH  EQUATION  1/ 

Donn  G.  DeCoursey£/ 

INTRODUCTION 

This  paper  presents  an  equation  that  defines  the  surface  runoff 
hydrograph  and  develops  a  method  for  determining  the  constants  in  the 
equation.   It  was  derived  in  studies  of  six  watersheds  in  the  Oklahoma 
portion  of  the  Washita  River  Basin.   The  equation  was  developed  as  one 
part  of  a  comprehensive  research  program  to  determine  the  changes  in 
flow  regimen  along  the  main  stem  of  the  Washita  River  associated  with 
the  application  of  flood  control  and  watershed  management  measures  on 
upstream  tributaries. 

It  is  believed  that  the  equation  fills  the  need  for  a  general 
functional  relation  that  can  be  used  in  the  electronic  computer  analysis 
of  streamflow  problems.  Many  conceptual  models  or  rational  approaches 
have  been  proposed  for  delineating  the  runoff  hydrograph:!.'  and  some  of 
these  concepts  have  been  used  in  developing  this  equation.   The  emphasis 
here,  however,  has  been  placed  on  having  the  equation  related  to  various 
hydrograph  properties. 


1/   Soil  and  Water  Conservation  Research  Division,  Agricultural  Research 
Service,  USDA,  in  cooperation  with  the  Oklahoma  Agricultural  Experiment 
Station. 

2^/  Hydraulic  engineer  (Research),  USDA,  Chickasha,  Okla. 

3/   Chow,  Ven  Te.,   1964.   Runoffo   Handbook  of  Applied  Hydrology. 
Section  14,  pp.  14-27  to  14-35.   McGraw  Hill  Book  Co.,  New  York. 


DERIVATION  OF  THE  EQUATION 

Perhaps  the  best  single  equation  for  representing  a  hydrograph  was 
derived  by  Edson^/  and  discussed  by  Bloomsburg^^/; 

q  =  a  tb  e-ct  [1] 

where  q  is  discharge  in  cubic  feet  per  second,  t  is  time  in  hours  from 
beginning  of  the  rise,  e  is  the  base  of  natural  logarithms,  and  a,  b,  and 
c  are  constants. 

The  equation  is  based  in  part  on  the  assumption  that  the  cumulation 
of  drainage  area  contributing  to  runoff  is  proportional  to  some  power  of 
time; 

A-^tb  b  >  1.0  [2] 

If  runoff  is  expressed  in  watershed  surface  inches,  q  will  also  be  pro- 
portional to  some  power  of  time; 

qc5<^tb  b  >  1.0  [3] 

A  third  assumption  is  that  the  reservoir  action  of  the  stream  causes 
the  discharge  to  decrease  exponentially: 

q^  e-ct  c  >  1.0  [4] 

Since  a  variable  which  is  proportional  to  two  functions  is  also  proportional 
to  their  product: 

qc><t^  e"'^*^ 

With  a^   as  the  proportionality  constant,  the  equation  becomes: 

q  =  a  t^  e-ct  [1] 


4/  Edson,  C.  G.   August  1951.   Parameters  for  relating  unit  hydrographs 
to  watershed  characteristics.   Trans.  Amer.  Geophys.  Union,   32(4):  591-596, 

5/   Bloomsburg,  G,  L.  March  30,  1960.   A  hydrograph  equation.   Agricul- 
tural Research  Service,  SWORD  Research  Report  No,  325,  13  pp. 


The  constants  in  this  equation  may  be  solved  for  a  particular 
watershed  in  terms  of  its  hydrograph  characteristics;  time  to  peak,  peak 
discharge,  and  the  volume  of  runoff.   This  was  done  for  the  hydrographs  of 
six  watersheds.   The  major  features  of  the  six  hydrographs  used  for  the 
analysis  are  presented  in  table  1—'.      Figure  12.'   shows  the  hydrograph  for 
Rock  Creek  and  the  hydrograph  as  obtained  from  the  equation. 

The  recession  side  of  the  equation  was  depleted  too  quickly.   This 
was  characteristic  of  all  six  watersheds.   Therefore,  it  was  decided  to 
combine  the  above  equation  with  a  second  equation  that  would  define  the 
slope  of  the  recession.   The  equation  used  for  the  recession  was: 

t-to 


q  =  Qo  e  [5] 

in  which  q  is  the  discharge  in  cubic  feet  per  second  at  time  t;  q^  is  the 
discharge  at  time  to;  e  is  the  base  of  natural  logarithms;  and  d  is  a  re- 
cession constant. 

Equation  5  may  be  rewritten  by  taking  the  natural  logarithm  of  both 


sides: 


In  q  =  In  qo  - 


t-to 


or 


1         ^  .  *^0     t 

In  q  =  In  qo  +  ^  -  ^ 

Since  qo  and  t^  are  constants  in  the  equation;  In  q  is  the  dependent 

o 
variable.  In  qg  +—  is  the  dependent  variable  intercept  for  t  equal  to 

zero,  t  is  the  independent  variable,  and  -  ■=.  is  the  slope  of  the  line.   For 

d 

a  given  watershed,  the  numerical  value  of  d  may  be  obtained  by  plotting  the 

normal  recession  for  the  surface  runoff  hydrograph  in  semilogarithmic  form 

and  drawing  a  straight  line  through  the  results.   The  slope  of  this  line 

1 
xs  --. 


6^/  All  tables  and  figures  are  presented  in  the  appendix, 


By  combining  equations  1  and  5  and  integrating,  the  total  volume  of 
runoff  during  the  period  of  the  hydrograph  is  given  by: 

Q  =  po  a  t^  e-c<^  dt  +  r°°  q^  e   ^   dt  [6] 

o  to 

Q  is  the  total  volume  of  runoff  in  the  hydrograph  expressed  in  cubic 
feet  per  second  hours;  t^  is  equal  to  kt  ,  the  product  of  a  factor  and  time 
to  peak;  qQ  is  the  discharge  in  cubic  feet  per  second  at  the  time  when  re- 
cession begins  (tg);  t  is  the  time  in  hours  from  the  beginning  of  the  rise; 
a,  b,  and  c  are  constants  obtained  by  solution  of  the  hydrograph  equation; 
and  d  and  k  are  constants  obtained  from  analysis  of  hydrograph  character- 
istics. 

The  first  part  of  equation  6  expresses  the  form  of  the  hydrograph 
through  the  peak  to  the  point  of  recession.   Therefore,  equation  1 

q  =  a  t^  e-ct 

may  be  used  to  help  solve  for  constants  a,  b,  and  c.   The  peak  flow  occurs 
where  the  first  derivative  of  q  with  respect  to  t  is  zero. 


^  =  a  tb  (-c)  e-ct  +  a  b  t^-l  e'Ct 
dt        ^      ^ 


and  by  simplifying, 


4a  =  _§_  (.tb  c  +  b  tb-1) 


dt    gCt 


— i   ia    pniifll     tn    zpro_       Since  — - 
c 
A 

tb  c  =  b  tb-1 


The  maximum  or  minimum  occurs  where  —^   is  equal  to  zero.   Since 

dt  gCt 

is  nonzero,  the  maximum  value  of  q  occurs  where: 


or  since  t  is  always  positive: 

tb    b 


tb-1  -  c 


=  -  =  t 


At  this  point  t  =  tp;  therefore: 


c  =  -^  [7] 


where  tp  is  the  time  to  peak  in  hours. 


By  substituting  t  =  tp  in  equation  1,  the  equation  for  peak  discharge 
is : 

qp  =  a  tpb  e-ctp 


letting  c  =  £_  and  solving  for  a: 
tp 


q£  e 


b 


a  =  -^i^  [8] 

t  " 

Thus,  the  constants  a  and  c  in  equation  1  can  be  expressed  in  terms  of  b, 
tp,  and  qp. 

Next,  q-  for  use  in  equation  5  is  determined  by  setting  t  =  tg  in 
equation  1: 

qo  =  a  t^b  e-cto  [9] 

By  substituting  Lq   =  ktp  and  equations  7  and  8  into  equation  9: 

b         .^^ 

qo  =  ^^   (kt_)b  e   tp 

t  b 
^P 


and  simplifies  to: 


q^  =  qp  k^  eb  e'^^  j^qj 


By  substituting  t^   =  kt-  and  equations  7,  8,  and  10  into  equation 
6,  the  total  volume  of  runoff  in  terms  of  b  is: 

b  bt 


Q  =  J^h  5£_!:  tb  e"  "^  dt 


b  tpb 


r  °°  „     ,,b  „b  .-bk 


t-ktp 


+   I        q^  k"^  e''  e-'^*^  e  d        dt  [n] 

'P 


J        > 
kt. 


Solution  of  equation  11  is  simplified  if  time  is  transformed  such 
that  tp  is  unity.   Therefore,  letting  t]^  =  — ,  when  t  =  t  ,  t]^  =  1.0, 
t  =  tj^tp,  and  dt  =  tp  dtj^;  equation  11  becomes: 

b  ^^l^P 


_  tit^-kt^ 

r    q  k^  e^  e-^k  e      d    t  dt, 
J     P  P   1 

kt-, 


By  simplifying  and  using  ti  =  1,0: 

P 


Q  =  qp  eb  tp  P^  t^b  e-bti  dt^ 


titp-ktp 
+  q  eb  kb  e-bk  t   p°e"     d    dtj^ 


By  dividing  both  sides  by  qpt   and  solving  the  integrals: 


ss 


.-bti   li_  . 


t^b   btjb-l   b(b-l)ti 


b-2 


-b    (-b)2       (-b)3 


(-b)b         (-b)b+l_ 


—  I  o 


t„k 


+  kb  e^  e-bk  e  d 


tit. 


~»  00 


or  by  simplifying,  the  second  half  of  the  equation  becomes: 


qptp 


.-bt- 


I-.  b 


-b 


bt 


b-1 


(-b)2 


b(b-l)t]^ 
(-b)3 


b-2 


«    •  •  9 


+  (-i)b-i^+  (-i)b_.±L^ 

(-b)b         (-b)b+l 


~ik 


—  I  o 


+  d_  kb  eb(l-k) 
t 
P 


[12] 


In  the  above  equation  all  things  are  known  except  b;  Q  is  the  volxmie 
of  runoff;  q^  is  the  instantaneous  peak  runoff  rate;  t  is  the  time  to  peak; 
d  is  the  average  recession  constant  for  the  watershed;  and  k  is  the  ratio  of 
the  time  at  which  the  recession  begins  to  the  time  to  peak. 

Since  there  is  no  direct  solution  of  equation  12  for  b,  a  graphical 
solution  is  necessary.   The  graphical  solution  would  be  a  family  of  curves 


in  which 


^■_  and  —  are  known  parameters.   Such  a  family  of  curves  would 
^P  P      P 


be  needed  for  each  value  of  k.   For  practical  purposes,  it  was  assumed  that 
one  family  of  curves  based  on  a  constant  value  of  k  would  be  satisfactory. 
In  order  to  establish  the  factor  k,  an  average  value  for  the  six  watersheds, 
about  1.8,  was  used  as  the  first  trial  value.   The  value  of  1.8  tended  to 
cause  the  recession  to  start  too  late  even  though  it  was  fairly  close  to  all 
actual  values.   This  may  be  partially  explained  by  the  fact  that  values  of 
the  equation  q  =  a  t"  e"^*^  are  falling  considerably  faster  on  the  recession 
side  of  the  hydrograph  than  are  the  values  of  the  actual  hydrograph.   For 


t-t, 


this  reason,  the  equation  q  =  q^  e 


becomes  effective  at  a  lower  value 


of  qo  than  it  should.   Further  analysis  indicated  that  this  is  partially 
corrected  by  assuming  that  the  recession  begins  at  an  earlier  period,  k  equal 
to  about  1.5.   The  value  of  1.5  appeared  to  be  satisfactory  for  all  six 
watersheds.   On  this  basis,  it  was  assumed  that  one  family  of  curves  based 
on  a  k  value  of  1.5  would  be  satisfactory  for  all  cases.   Figure  2  shows  the 
hydrograph  of  Rock  Creek  and  two  equations  obtained  by  using  values  of  1.5 
and  1.8  for  k. 


By  substituting  the  value  1„5  for  k  in  equation  12,  the  equation  is: 


-±-=   eb 

q  t 
p  p 


>-bti  |"ll^  .  ^<^l^'^  +  b(b-l)tib-2 

f-^   (-b)2    (.b)3 


1.5 


+  (.l)b-l  il£l  +  (.l)b  _bl      j      +1-  (i.5)b  e-0.5  b 


.b    ^  ^   /_KNb+l 


(-b)°         (-b)"^^   I        ^p 


f-   (1.5)°  e-^-^  °     [13] 


— '  o 


The  family  of  curves  for  solution  of  equation  13  is  shown  in  figure  3, 

— ~—  was  plotted  versus  b  for  various  values  of  — , 
^P^P  *^p 

For  constant  values  of  k,  equation  13  may  be  shown  as: 

-^=  f(b)  +±-   g(b)  [14] 

qptp  tp 

In  case  the  ranges  of  ^   ,  — ,  or  b,  shown  in  figure  3,  are  not  wide 

qptp  tp 

enough,  values  of  f(b)  and  g(b)  for  k  =  1.5  and  b  =  0  to  40  are  presented  in 
table  2,   These  values  may  be  used  to  calculate  the  points  required  to  plot 
a  family  of  curves  similar  to  that  of  figure  3.   They  may  also  be  used  in 
graphical  interpolation  for  precise  values  of  b  by  expanding  the  graph  of 

— -* —  versus  b  for  a  particular  value  of  — ,   The  set  of  four  integer,  values 
qptp  tp 

of  b  are  estimated  from  the  family  of  curves  such  that  they  bracket  the  true 

value  of  b.   Using  these  four  values  and  — ,  values  of  ^^  are  calculated 

t  at 

P  ^P  P 

and  plotted  versus  b  on  a  large  scale.   The  points  will  plot  as  a  smooth 

curve  (almost  linear).   The  exact  value  of  — ^ —  is  entered  and  the  value  of 

Q  t 
^P  P 

b  is  found o 


10 


Use  of  the  hydrograph  equation  begins  with  solution  of  the  variable 
b.   Q,  Qn,  t  ,  and  d  are  determined  for  a  hydrograph;   b  is  then  obtained 
from  the  curves  of  figure  3,  or  perhaps  more  accurately  by  using  graphical 
interpolation  with  the  data  presented  in  table  2.   Variables  a  and  c  are 
calculated  from  equations  8  and  7,  respectively.  With  all  the  variables 
known,  the  discharge  at  any  time  can  be  calculated  from  equations  1  and  5. 
Equation  1  defines  the  hydrograph  from  time  t  equal  to  zero  to  tp  (1.5  times 
the  time  to  peak).   Equation  5  defines  the  recession  of  the  hydrograph  from 
time  to  to  the  end  of  the  hydrograph.   Before  equation  5  may  be  used,  the 
value  of  qo  at  time  t^   must  be  calculated  from  equation  1.   Table  3  is  a 
sample  calculation  of  the  hydrograph  of  Rock  Creek.   Figure  2  is  a  plot  of 
the  hydrograph  and  the  equation  for  the  Rock  Creek  watershed. 

Figures  4  through  8  are  the  hydrographs  and  equations  for  the  other 
five  watersheds. 
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APPENDIX 


TABLE  1. "-Major  features  of  hydrographs  for  six  tributaries 
of  the  Washita  River,  Okla. 


:  Time; 

Drainage  : 

:   to  : 

Peak  : 

Rec,   : 

Time 

Name      ; 

area   : 

Volume   : 

peak: 

rate  : 

constant: 

constant 

•  Sq.  miles: 

C.f.s.  hr. 

•  Hr.: 

C.f.s.: 

d    : 

k 

Barnitz  Creek 

243   : 

156,815 

>   21 

>   9,500" 

4.82  \ 

1,52 

Little  Washita 

227   : 

146,490 

:   14 

:  14,000: 

6.11  : 

•   1.57 

Pond  Creek 

319   : 

205,860 

:   18 

:  18,000 

6.55 

1.61 

Rock  Creek 

:    138 

89,056 

I       7 

:  13,000 

:   5.89 

:   1.86 

Rush  Creek 

\         206 

•   132,938 

:   4 

:  20,000 

:   4.67 

:  2,00 

Sugar  Creek 

:    211 

\      136,158 

:  10 

:  11,000 

:   5.39 

:   1.90 
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TABLE  2. — Values  of  f (b)  and  g(b)  for  use  in  equation  14 


b 

;     f(b) 

:    g(b) 

;    b 

:     f(b) 

:     g(b) 

0 

:   1,50000 

:    1,00000 

:   21 

:   0.53188 

:   0.13735 

1 

:   1.20196 

:   0.90980 

:   22 

:    .52137 

. 12496 

2 

i   1.06552 

.82773 

;   23 

.51140   ! 

.11369 

3 

:   0.97854 

.75306   : 

:   24   i 

:    .50192   ! 

.10343 

4 

:    .91487   : 

.68513   i 

•   25 

:    .49290   : 

.09410 

5 

\          . 86462   : 

.62333   : 

26   : 

.48430   ! 

.08561 

6 

.82306   : 

.56711   : 

27   : 

.47610   : 

.07789 

7   : 

.78761   J 

.51595   : 

28   i 

.46827   : 

.07087 

8   : 

.75667   : 

.46941   \ 

29   \ 

.46077   ; 

.04447 

9   \ 

.72924   \ 

.42707   : 

30   : 

:    .45360   : 

.05866 

10   : 

.70461   : 

. 38854   : 

31    ! 

.44672   i 

.05337 

11   : 

.68227   : 

.35350   : 

32   i 

.44012   i 

.04855 

12   : 

.66185   : 

.32161   : 

33   : 

.43378   : 

.04417 

13   : 

.64306   : 

.29260   ! 

34   i 

.42769   : 

.04019 

14   i 

.62569   i 

.26620   : 

35   : 

.42183   i 

.03656 

15   : 

.60955   i 

.24219   ■ 

36   i 

t    .41619   : 

.03326 

16   i 

.59449   : 

.22035   ; 

37   i 

.41076   : 

.03026 

17   ': 

. 58040   : 

.20047   : 

38   i 

.40552   i 

.02753 

18   J 

.56717   i 

. 18239   J 

39   i 

.40046   ! 

.02505 

19   : 

.55472   i 

.16593   : 

40   i 

:    .39558   i 

.02279 

20   i 

,54298   ; 

.15097   : 
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TABLE  3. — Hydrograph  calculation  for  Rock  Creek 


Equations  used:   q  =  a  t^  e~^*^ 
Constants:        Q  =  89,056 


qp  =  13,000 


t-kt, 


q  =  qo  e 


d  =  5.89 
k  =  1.5 


S 


=  7 


ktp  =  10.5 


q^tp    (7) (13,000) 


?i^056 ^  0^9786 


d_^  5^^  0.8414 
t     7 
P 


b  =  11.00  from  figure  2 

=  11.027  by  graphical  interpolation 

c  =  b_=  iJU027  ^  1^57529 


a  = 


-  ^P  ^^  _  (13.000)  ell.027 


7II.O27 


log  a  =  log  13,000  +  11.027  log  e  -  11.027  log  7 
=  -  0.415993  =  9.584007  -  10.0 
=  0.383713 


q  =  a  t^  e-ct 
log  q  =  log  a  +  b  log  t  -  ct  log  e 

=  -  0.415993  +  11.027  log  t  -  0.684139t 


1 
2 
3 
4 
5 
6 


lofi  t          ! 

!      log  q        : 

1 

)           : 

:    -1.100132 

0 

.301031     ! 

!     1.535187 

34 

.477121      i 

:     2.792803 

621 

.602060      ! 

'.            3.486367 

3,065 

.698970      ! 

'     3.870854      i 

7,428 

.778151      i 

:     4.059844      ; 

11,477 
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Table  3. --Continued 


t 

:      lofi  t       : 

log  q 

1 

7       : 

i      .845098      ! 

4.113930 

:   13,000 

8       : 

:      .903090      j 

4.069268 

11,730 

9       : 

.954243      ! 

3.949194 

:    8,896 

10       : 

1.000000      : 

3.769617 

5,883 

10.5     : 

1.021189      I 

3.661199      ! 

4,584 

t-ktp 

q  =  <lo  «    ^ 

log  q  =  log  qo  - 


t-kt^ 


log  e 


=  3.661199  -  ^"^Q-^  (0.434294) 
5. 89 


=  4.435407  - 

-  0.073734t 

11       : 

:     3.624333 

:    4,211 

12       i 

!     3.550599 

i    3,553 

14       i 

:     3.403131 

\          2,530 

16       i               i 

\            3.255663 

:    1,802 

18       i               i 

:     3.108195 

1,283 

20       i 

:     2.960727 

:      914 

24       i 

:     2,665791 

I      463 

28       i 

:     2.370855 

235 

34       i               i 

!     1.928451 

85 

40       i               i 

:     1.486047      • 

I       31 

50       :               ! 

!      .748707      '. 

6 

15 


14  r 


(O 


o 

o 
o 
o 


u 
o 
a: 

< 

X 


24 

TIME   IN  HOURS 


36 


42 


Figure  1.— Rock  Creek  hydrofraph;  equation  q 
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Figure  2.—  Rock  Creek  hyixmgrmfh. 


17 


pnTTn^iMiiii|iiiiinTT|TTTTTTTTT|TTTTTTOTTTTTT7TTTJTWTTTTT]TTTTH'lll|lMllllll|l|||||||^ 


iiiiiiirliinrrniliiiiiiiiiliiiiiiiiiliiiiiHiil Iiiiiiiiiil iiiiliiiiiiiiiliiiiiiiiiliiiiiiiiiliiiiMiiiliiiiiiiiiliiiiliiiiliiiiiiiii 


^ 


in 


to 


o 


lO 


d.d 


i'b. 


in 
6 


18 


14 


12 


10 


EQUATION 


HYDR06RAPH 


18  24  30 

TIME    IN  HOURS 


48 


Flfur*  4.»Bamlts  Creek  hydrograph. 
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risur«  3.— Llttl*  HMhlta  RItmt  hydrograph. 
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Figure  6.— Pond  Creek  hydrogreph. 
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Flgur«  7.— Rush  Creek  l|ydrofr*ph. 
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Figure  8.— Sugar  Cre«k  hydrogr«ph. 
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